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ABSTRACT 

We present evidence that the sound-speed variation with solar activity has a two-layer config- 
uration, similar to the one observed below an active region, which consists of a negative layer near 
the solar surface and a positive one in the layer immediately below the first one. Frequency dif- 
ferences between the activity minimum and maximum of solar cycle 23, obtained applying global 
helioseismology to the Michelson Doppler Imager (MDI) on board SOHO, is used to determine 
the sound-speed variation from below the base of the convection zone to a few Mm below the 
solar surface. We find that the sound speed at solar maximum is smaller than at solar minimum 
at the limit of our determination (5.5 Mm). The min-to-max difference decreases in absolute 
values until ~ 7 Mm. At larger depths, the sound speed at solar maximum is larger than at solar 
minimum and their difference increases with depth until ~ 10 Mm. At this depth, the relative 
difference (Jc^/c^) is less than half of the value observed at the lowest depth determination. At 
deeper layers, it slowly decreases with depth until there is no difference between maximum and 
minimum activity. 
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The Sun has an approximately 22-year mag- 
netic cycle, where the dipolar magnetic field at 
the solar poles reverses each 11 years or so. Dur- 
ing the reversal, the magnetic activity is at a mini- 
mum and very few sunspots are visible on the Sun, 
sometimes none can be seen. At the maximum in 
the cycle, the number of sunspots visible on the 
Sun can be more than 100 at one time. A typical 
sunspot has a lifetime of a few weeks. The radio, 
ultraviolet and X-ray emission of the Sun also in- 
crease significantly during solar maximum. The 
total solar irradiance too has a strong correlation 
with the sunspot number. Significant flares occur 
at all phases of the sunspot cycle, but the num- 
ber of M-class and X-class fiares tends to follow 
the sunspot number. There are several proxies 
of solar activity currently used, probing the so- 
lar atmosphere at di f ferent heights (se e, for exam- 
ple, I Jain et allliooi iHathawavllioiol ). The most 
accepted scenario is that the solar cycle is pro- 



duced by dynamo processes within the Sun and it 
is driven by differential rotation and convection. 

As initially reported by IWoodard fc Nove^ 
()l985l) . it is now well established that the fre- 
quency of the solar acoustic modes is strong 
correlated with solar activity. This has been 
observed for modes with degree up to 3 (e.g., 
Salabert et all |2004 'Chaplin et al' '2007|) up 



to about 200 (e.g., .gachma nn fc Brown , 1993t 



Jain fc Bhatnagai|2003l:lAntiall2003HDziembowski fc Goode 



20051: IJain et all ^2009^ and up to 100 
Rhodes et al.l 2002; Rabello-Soares et a 



Rabello-SoaresI 1201 it iRhodes et all l201l[ ). The 



2om 



observed mode frequencies can be inverted to ob- 
tain the sound speed in the solar interior where 
the waves travel through. However, in analyzing 
the changes in mode frequency with the mag- 
netic solar cycle, we have to keep in mind that 
magnetic fields can change not only the thermo- 
dynamic structure of the medium the waves travel 
through, which, in turn, changes the frequencies 
of the waves but also the plasma waves can be 
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directly affected by the magnetic fields through 
the Lorentz force. In short, the modification to 
the acoustic mode frequencies results from both 
structural and non-structural e ffects of the mag- 
netic fields (see Lin et al. 2006[) . 

The measured solar-cycle frequency shifts de- 
pend strongly only on the frequency of the mode 
afte r they have be en weig hted b y the mode iner- 
tia ( Libbrecht &: Woodard 1990l) . This indicates 
that the dominant structural changes during the 
solar cycle, so far as they affect the mode frequen- 
cies, occur near the surface. At the outermost 
layers, the timescale of the heat transfer is smaller 
than the period of the oscillation and they are 
not adiabatic anymore. At the moment, there 
is no general agreement as to the precise phys- 
ical mechanism near the surface that gives rise 
to the freque ncy variation, (see, for example, 
Li et al. Several authors confirmed that all 

or most of the physical changes were confined to 
the shallow layers of th e Sun (e.g., [Basu 200 j: 
Eff-Darwich et all l2002t iDziembowskTfc GoodA 



20051 : IRabell o-Soares & KorzennikI l2009l etc). 
However, iBaldncr & Basu (2008) observed a 
small, but significant change in the sound speed 
at the base of the convective zone and for r > 0.9 
Rq. a ch ange in the inte ri or str ucture was also 
found by iBasu fc Mandell (|2004l ). specifically a 
variation in adiabatic index, Fi, near the sec- 
ond helium ionization zone (~ 0.98 Rq). On 
the other hand, changes in zonal and meridional 
flows correlated with solar activity have been ob- 



serve d by several authors (jScb ou 1999: Howe et al 



Basu &: Ant ial l2006l etc) 



2OOOI: IVorontson et all l20CI2t iHowe et all l2005t 



Local helioseismology techniques have shown 
that below an active region there is a two-layer 
structure with a negative variation of the sound 
speed in a shallow subsurface layer and a positive 
variation in th e deeper interior relative to a quiet 
Sun region (see Kpsovichev et al.'2000'; Basu et aL 
2004 iBogart etal., ,2008.: Kosovichev et al.ll2011 , 
and references within). A similar variation has 
been observed for the adiabatic in d ex, Fi below 
an a ctive region ( Basu et al. 20041 : Bogart et al. 
20081 ). Because global activity levels are much 



smaller than in active regions, a similar structure 
variation would be much smaller in a global scale, 
if it occurs at all. 

Rabello-Soared ( 2011 ') estimated the frequency 



differences between the activity minimum and 
maximum by carefully fitting the variation of 
the acoustic mode frequency with solar activ- 
ity for modes with 20 < ^ < 900 observed by 
MDI/SOHO. Here, these frequency differences 
(described in Section 2) are inverted to look for 
the interior sound-speed variation with solar cycle. 

2. Data used 

Full-disk Doppler images obtained at a one- 
minute cadence by MDI Dynamics and Structure 
observing modes were analyzed (I Scherrer et al 



Il995[ ). The first one has higher spatial resolution 
(1024^ pixels) and is available every year for two 
or three months of continuous data. Data for the 
years 1999 to 2008 were used. The length of the 
time series varies from 38 days (in 2003) to 90 days 
(in 2001). The second observing mode is available 
year around, but is subsampled in order to fit the 
limited telemetry (modes up to Z w 250). It is di- 
vided in 72 day time series and data from early 
1996 to April 2010 were used. They will be called 
the Dynamics and Structure sets from now on. 
Each velocity image is decomposed into spherical- 
harmonic components (iSchou et al...2002, ). 

Two different peak-fitting algorithms were used 
to fit the power spectra and obtain the mode fre- 
quencies. One of them known as the MDI peak- 
fitting method is described in detail by ISchou 
(1992) (see also lSchou et al.ll200"ilLarson fc Schon 
I2OO8,, ,2009! ). A spherical-harmonic decomposition 
is not orthonormal over the solar surface that can 
be observed from a single view point resulting in 
what is referred as spatial leakage. At high de- 
grees, the spatial leaks lie closer in frequency (due 
to a smaller mode separation) resulting in the 
overlap of the target mode with the spatial leaks 
that merges individual peaks into ridges, making 
it more difficult to estimate unbiased mode fre- 
quencies. The MDI peak-fitting method works 
very well for low and medium-Z modes where the 
leaks are well separated from the target mode. 
For high-Z modes, which have spatial leaks over- 
lapping t he target inode and blended in a sin- 
gle ridge, iKorzcnnik' (199 81) method wa s used (for 
details, see iRabello-Soares et al.l 120081 ). The fit- 
ting was carried out only for every tenth I for 
modes with 100 < I < 900. This method was ap- 
plied only to the Dynamics time series which were 
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Fourier transformed in small segments (4096 min- 
utes) whose spectra were averaged to produce an 
averaged power spectrum with a low but adequate 
frequency resolution to fit the ridge while reduc- 
ing the realization noise. These two peak-fitting 
methods will be called ML and HL respectively. 
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Fig. 1. — Mean radio flux for the observational 
periods in the Dynamics set (in black) and in the 
Structure set (in red). The error bars correspond 
to the error of the mean. 

The ML method was applied to both observ- 
ing sets and the HL only to the Dynamics set. In 
summary, three types of data were used in this 
work: (1) frequencies obtained by both ML and 
HL methods applied to the Dynamics set ("Dy- 
namics ML-fHL"); (2) frequencies obtained by 
only the ML method applied to the Dynamics set 
("Dynamics ML"); and (3) frequencies obtained 
by only the ML method applied to the Structure 
set ("Structure ML"). 

The estimated mode frequencies, Vnim, were 
parameter ized in terms of C l ebsch -Gordan co- 
efficients (jRitzwoller fc Lavelvl Il991l ) . The cen- 
tral frequency, i.e., the frequency free of split- 
ting effects, is taken to be the frequency given 
by m = in the parameterization. For each of 
these three sets, the central frequency, Vni, for 
each time period was fitted assuming a linear and 
a quadratic relationship with the correspondent 
solar-activity index and using a weighted least- 
sq uares minimization. F or a detailed description, 
see lR,abello-Soar"e3(|2nilh . The solar radio 10.7 cm 
daily flux (NGDC/NOAA) was used as the solar- 
activity proxy (Figure [l} . 

Although, there is a very high linear correla- 



example. IChaplin et al.l 120071 ). More recently. 



Rabello-Soared (|201ll ) found some evidence of a 
quadratic relationship indicating a saturation at 
high solar activity (^ 350 x 10^^^ J s^^ m~^ 
Hz~^). A similar effect has been seen in frequen- 
cies at activity regions with a large surface mag- 
netic field using ring analysis (Basu et al. 2004). 
The frequency shifts obtained using the quadratic 
fitting were used here. 

Modes with a negative linear coefficient in the 
quadratic fitting were excluded from the present 
analysis. They also have a positive quadratic co- 
efficient, indicating that they are less affected by 
low and medium levels of solar activity. They 
correspond to 5%, 11% and 16% of the modes 
of the Structure ML, Dynamics ML and Dynam- 
ics ML-I-HL sets, respectively. In the Dynamics 
ML-I-HL mode set, the excluded modes make for 
the lack of modes in the middle of the ^ — di- 
agram in Figure [2l The /-modes obtained with 
the HL method (about 20 modes) have a distinc- 
tive behavior from thos e obtained with the ML 
method (see Figure 7 in iRabello- Soared 2011 ) . It 
is not clear if this behavior is an artifact of the fre- 
quency determination or not. They were excluded 
from our analysis. 
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Fig. 2. — Modes used in the inversion for the Dy- 
namics ML-f-HL, Dynamics ML and Structure ML 
sets are shown as black squares, red crosses and 
blue circles, respectively. 

The solar activity minimum-to-maximum fre- 
quency shift 5v^^i is defined as the difference be- 
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tween the fitted frequency at the maximum and 
at the minimum of solar activity. The minimum 
and maximum activity correspond to the Dynam- 
ics observing periods centered around April 2008 
and May 2002 respectively (Figure [J). These ob- 
serving periods were also used to calculate the fre- 
quency difference for the Structure set. The cor- 
responding radio fluxes are 71 and 185 J s~^ m~^ 



Hz 



for both sets. The minimum-to-maximum 



frequency shifts Sv^'f obtained using a linear or a 
quadratic fit to the solar activity are very similar. 

3. The surface term 

The scaled frequency differences Qni 5vni / Vni 
(between two models or between the Sun and a 
model) can be expressed as Hi{vni/L) -\- H2{vni), 
where Hi is the contribution of the interior sound- 
speed difference, H2 of the differences in the sur- 
face layers, Qni is the mode inertia normalized by 
the inertia of a radia l mode of the same frequency 
and L — -\- 1) ( Christensen-Dalsgaard et al.l 
[1989,) ■ H2 is also called the surface term, Fgurf ■ 
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Fig. 3. — Scaled minimum-to-maximum fre- 
quency differences obtained using Dynamics 
ML-I-HL (black squares), Dynamics ML (black 
crosses) and Structure ML (small black circles) 
sets. Green, red and blue circles correspond to 
the averaged values for Dynamics ML-f-HL, Dy- 
namics ML and Structure ML, respectively. The 
solid lines are the interpolated splines representing 
the surface term, -Fsurf- 

Figure [3] shows the scaled frequency differ- 
ences between solar minimum and maximum, 
Qni^i^^i/^nii obtained using Dynamics ML+HL 
(black squares). Dynamics ML (black crosses) and 



Structure ML (small black circles) sets. They are 
basically a function of frequency. The near-surface 
contribution can be removed by subtracting a 
smooth function of frequency, i^surf(i^), from the 
scaled frequency shifts, thus allowing the search 
for structure variations in the solar interior with 
the solar cycle. Qni Sv^l/vni were averaged over 
200-/iHz intervals. In each interval, outliers large 
than 3-cr were removed and the weighted average 
calculated. The results are represented in Fig- 
ure [3] as green, red and blue circles respectively 
for Dynamics ML-I-HL, Dynamics ML and Struc- 
ture ML sets. A cubic spline interpolation was 
applied to the 17 and 21 frequency- knots for the 
ML and ML-I-HL sets, respectively (solid lines 
in Figure [S]) . Figure H] shows the residuals after 
subtracting the surface term given by the splines. 
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Fig. 4. — Residuals, QniSv^i jv^i, obtained after 
subtracting Fsurf. The black squares, red crosses 
and blue circles correspond to Dynamics ML-I-HL, 
Dynamics ML and Structure ML, respectively. 



4. Results 

A nonasymptotic inversion technique can be ap- 
plied to infer the sound speed (and density) in- 
side the Sun using the frequency differences be- 
tween the Sun and a reference solar model (e.g. 
Antia fc BasijlTggi . and references within). Here 
the solar-cycle frequency shifts, Sv^^^ , were in- 
verted to determine the differences in sound speed 
between minimum and maximum solar activity in- 
side the Sun. The inversion has been carried out 
using the multi plicative optim ally localized aver- 
ages technique (jGough Il985h . The relative fre- 
quency shifts bv^^ jvni are plotted in Figure [S] 
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Fig. 5. — Relative frequency differences between 
solar maximum and minimum, after removing the 
surface term, used in the sound-speed inversion. 

The inferred sound-speed differences are shown 
in Figure |6l The results obtained using Dynamics 
ML-I-HL (black). Dynamics ML (red) and Struc- 
ture ML (blue) sets follow the same general trend. 
The sound-speed difference increases slowly with 
radius until about 0.98 Rq after which it decreases 
sharply. The sound-speed difference is larger than 
6-cr at its maximum for all three sets. The inver- 
sion results are identical if a linear fit of the mode 
frequencies to the solar-activity index is used in- 
stead of a second-degree polynomial. 

The solution errors for the Dynamics ML+HL 
set are smaller by as much as 25% than for the 
other two sets for r > 0.85 Rq and 15% for 
r < 0.85 Rq. Also, the averaging kernels for Dy- 
namics ML+HL are slightly more localized close 
to the surface than the other two sets due to 
the inclusion of high-^ modes. In Figure |6l the 
averaging-kernel widths (defined as the difference 
between the first and third quartile points) are 
given by the horizontal error bars for a few se- 
lected locations. In the top panel, the horizontal 
bars are the same for all three sets. In the bottom 
panel, they are in black for Dynamics ML+HL and 
in red for Dynamics ML and Structure ML sets. 



Fig. 6. — Inferred sound-speed differences ob- 
tained using Dynamics ML+HL (black circles). 
Dynamics ML (red line) and Structure ML (blue 
line) sets. The bottom panel shows the solu- 
tion near the surface in more detail. The vertical 
dotted lines mark approximately the regions that 
correspond to positive, zero and negative sound- 
speed variations below an active region in rela- 
tion to a quiet reg ion obtained by ring analysis 
(iBogart et al.ll2008l ). 



since they are very similar. Nevertheless, the av- 
eraging kernels are well localized and have almost 
no side lobes for the inferred solution from 0.53 to 
0.992 Rq for all three sets (Figure [71). The cross- 
term kernel measures the influence of the contri- 
bution from density on the inferred sound speed 
and it is insignificant in this range (bottom panel) . 
We did not include modes with Z < 20 in our anal- 
ysis since we are not interested here in the solar 
deeper layers. For r > 0.992 Rq, the averaging 
kernels become increasingly asymmetric and de- 
velop large side lobes. 

To check that the surface term is properly re- 
moved and is not causing the observed sound- 
speed variation, the surface term was estimated 
using four times more frequency-knots than used 
in Figure |3] (i.e., using 50 /iHz intervals). In this 
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Fig. 7. — Top: Samples of averaging kernels at 
five diff'erent radii near the surface for the Dy- 
namics ML+HL set. The target radius is shown 
by a vertical dotted line. Bottom: Averaging ker- 
nels at r = 0.992 Rq for Dynamics ML-I-HL (solid 
line), Dynamics ML (dashed line) and Structure 
ML (dotted line) sets. The cross term kernel is 
also plotted for all three sets. It is multiplied by 
100 and displaced vertically by -30 for visualiza- 
tion purpose. 



case, i^surf is not a smooth function of frequency 
anymore. The inversion of the residuals is still 
similar to Figure |6l The maximum sound-speed 
variation around 0.98 Rq is still statistically sig- 
nificant, it is larger than A-a (larger than 6-cr for 
the Dynamics ML-I-HL set). 

Figure [5] shows the sound-speed difference be- 
tween two different solar models (solid red line) 
and the correspondent inversion results (black cir- 
cles) obtained using the same mode set as our ob- 
servations given by Dynamics ML+HL. The inver- 
sion results are very similar using any of the other 
two sets. The two solar models are Model S and 
a model with identical physical assumptions but 
with a lower solar age, 4.52 Gyr, inste ad of 4.6 
Gyr ()Chr istcnscn -Dalsgaard et al.lll996 ). The in- 
version results agree very well with the true model 
difference for r < 0.992 Rq, specially taking into 
account the averaging kernel (dashed red line in 
the bottom panel). 
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Fig. 8. — Sound-speed differences between two so- 
lar models through most of the solar interior (top) 
and, in detail, near the surface (bottom). The 
solid red line is the true difference between the 
models, while the dashed red line (bottom panel) 
is the difference weighted by the averaging kernels. 
The black circles represent the inversion results us- 
ing the observed mode set for Dynamics ML+HL 
and the solid black line the inversion results with 
random noise added to the theoretical frequencies. 
The two dotted black lines correspond to the av- 
erage of 1000 realizations plus and minus 2.5-cr. 



Howe &: Thompson ( 19961 ) showed that random 



noise in the observed mode frequency will in- 
troduce spurious oscillatory structures in the so- 
lution as a result of the correlation of the er- 
rors in the solution at different radii. The in- 
version result at two different radii makes use 
of similar data sets and their correspondent er- 
rors. The solution errors are correlated even be- 
tween points farther apart than the averaging- 
kernel width. F or examples of sound - speed error 
correlation, see Rabello-Soares et al. ( 19991 ) and 
•B asu et all (|2003l ). To illustrate the effect of the 
error correlation in the inversion results presented 
here, it was carried out inversions of the theoreti- 
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cal frequency differences between the two models 
with noise added to them. To each frequency dif- 
ference the noise was generated as normally dis- 
tributed random numbers with the same standard 
deviation as the error in the observed data. The 
solid black line in Figure [8] shows the result for one 
realization. The spurious oscillatory structures in- 
troduced by the correlated errors, also present in 
Figure [6l are clearly seen. The two dotted lines 
correspond to the average of 1000 realizations plus 
and minus 2.5-(7, showing that the overall results 
are consistent with the estimated formal errors. 

To remove the spurious oscillation in our sound- 
speed determination, a sum of two Gaussian func- 
tions was arbitrarily chosen as a smooth function 
of radius, 5c^ /c^{r). It was 'weighted' by the av- 
eraging kernels and fitted to the inversion results 
of the two Dynamics sets using a nonlinear least- 
square fit (red solid line in Figure [9]). According 
to this simple model, the sound-speed difference 
increases with radius from around the base of the 
convection zone until 0.985 Rq (i.e., 10 Mm below 
the surface) at the second helium ionization zone 
where 8(? jc? = 3.9 x 10^^. Closer to the surface, 
it decreases steeply with radius. It is zero at 0.990 
i?0 (i.e., at a depth of 7 Mm) and, at 0.992 Rq 
(or 5.5 Mm deep), 6(? jc? — -9.2 xlO^^, where the 
sound speed at solar maximum is smaller than at 
solar minimum. 

5. Discussion 



Usi ng global helioseismology, iBaldner fc Basu 
POOSI ) found an increase with radius in the sound- 
speed variation with solar cycle for 0.86 < r/i?0 < 
0.95 (at the equator and at 15° in latitude) that 
agrees with our results within less than \.h-o. 
Here, we found that the variation in sound speed 
increases until 0.985 i?© after which it decreases 
fast with radius becoming negative for r > 0.990 
Rq. Posing the question of what is the physical 
mechanism behind this suddenly change in behav- 
ior at this particular depth. 

This crossing of the sound speed at solar max- 
imum from larger to smaller than at solar mini- 
mum is consistent in location with the crossing of 
the sound speed below an active region from larger 
to smaller than at a quiet region obtained by ring- 
diagr am analysis (jBasu et al. 1 12004 iBogart et al.l 
20081 ). The vertical lines in Figures [6] and [9] mark 
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Fig. 9. — Black circles and diamonds are the in- 
version results using Dynamics ML-I-HL and Dy- 
namics ML respectively. They are identical to the 
black and red circles in Figure [H The red line is 
the smooth function, Sc^/c^, which was weighted 
by the averaging kernels and fitted to the inversion 
results of the two Dynamics sets. The blue line is 
8<? I(? weighted by the averaging kernels. 



approxi mately the regions t hat correspond to Fig- 
ure 9 in lBogart et al.l (|2008l ) for positive, zero and 
negative sound-speed variations below an active 
region in relation to a quiet region. This agree- 
ment indicates that the sound-speed variation in 
the solar interior with solar cycle might be an over- 
all effect of the active regions local perturbation, 
and not a global structure change with solar activ- 
ity. The estimated sound-speed variation obtained 
here is based on averages over the whole Sun at a 
particular depth and over one or more solar rota- 
tions. 



The number of active regions is very well cor- 
related with the solar activity cycle. Using the 
data provided by SIDC (Solar Influences Data 
Analysis Center of the Royal Observatory of Bel- 
gium), the daily sunspot number averaged over 
the maximum-activity observing period (centered 
around May 2002) is about 120, while during the 



7 



minimum (April 2008) is less than 10. According 
to ring-diagram analysis, the sound-speed varia- 
tion at a depth of 14 Mm (i.e., 0.98 Rq) below 
a 16° X 16° region on the solar surface (tracked 
for 5.7 days) with a strong active region and with 
a magnetic activity index (MAI) of about 80 G 
i n comparison to a quiet region is: fx? jc? ~ 0.01 
( Basu et al. 20041 ). For a smaller sunspot with an 
MAI of 40 G, it is 0.05. In brief, the MAI is the 
average of MDI magnetograms over the same area 
in the solar surface and observing tim e as the an- 



alyze d region (for a full definition see iBasu et al 



20041 ). Performing a back of the envelope calcula- 



tion, if there are 11 regions with sizes of 16° x 16° 
(in the front and back of the Sun) with an MAI 
> 80 during solar maximum (or a larger number of 
weaker regions), the mean sound-speed variation 
will be 8(? jf? = 4 X 10""*, i.e., the value estimated 
here at this depth for the maximum-to-minimum 
solar-cycle variation. 

Time distance is another local helioseismol- 
ogy technique that provides information about 
the subsurface structure of active regions. While 
ring-diagram technique is based on the analysis 
of three-dimensional power spectra in a small re- 
gion of t he Sun to measure the oscillation fre- 
quencies (jHilll 119881 ). as in global helioseismology, 
time distance measures the wave travel times be- 



tween different points on the surface (jPuval et al 
[1991) • 

Both the ring-diagram and time-distance 
techniques give qualitatively similar results, re- 
vealing a characteristic two-layer seismic sound- 
speed structure. However, the transition between 
the negative and positive variations occurs at dif- 
ferent depths: ^ 4 Mm for the time-distance re- 
sult, and ^5 — 8 Mm for the ring-diagram inver- 
sions. It seems that this can be explained, at least 
in part, by differences in the se nsitivity and res- 



olutio n of the two methods (see iKosovichev et al 
201l[ ) 



Recently, lllonidis et ahl (|201ll) using time- 
distance technique detected strong acoustic travel- 
time anomalies around 65 Mm (i.e., at 0.91 Rq) 
below the solar surface where a sunspot emerged 
1 to 2 days later. These travel-time perturbations 
are an indication o£ sound-speed variations. In 
our simple model, 5?/? = 2.1 x 10~^ at this 
depth. However, the observed signal lasts only 
for ~ 8 hours and is not as strong in the layers 
immediately above the detected signal. 



In conclusion, using acoustic mode frequencies 
observed by MDI/SOHO through most of solar cy- 
cle 23, we have estimated the sound-speed differ- 
ence between solar maximum and minimum from 
0.53 to 0.992 Rq. The sound-speed difference 
seems to increase with radius from the base of the 
convection zone or above until 0.985 Rq (i.e., 10 
Mm below the surface) at the second helium ion- 
ization zone. Closer to the surface, it decreases 
steeply with radius. It is zero at 0.990 Rq (i.e., at 
a depth of 7 Mm) and, at 0.992 Rq (or 5.5 Mm 
deep), '5?/? = -9.2 x 10^''. The transition of 
the sound speed at solar maximum from larger to 
smaller than at solar minimum takes place at a 
similar depth as the transition of the sound speed 
below an active region from larger to smaller than 
at a quiet region obtained by ring-diagram anal- 
ysis, suggesting that the sound-speed variation in 
the solar interior with solar cycle could be an over- 
all effect of the active regions local perturbation, 
and not a global structure change with solar ac- 
tivity. 

We thank Jesper Schou and Tim Larson for 
providing the spherical-harmonic decomposition 
of the MDI images and the frequency tables from 
the MDI peak-fitting method. SOHO is a project 
of international cooperation between ESA and 
NASA. 
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